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We use a many-body rate-equation approach to calculate the thermopower of a quantum dot in 
the presence of an exchange interaction. At temperatures much smaller than the single-particle level 
spacing, the known quantum jumps (discontinuities) in the thermopower are split by the exchange 
interaction. The origin and nature of the splitting are elucidated with a simple physical argument 
based on the nature of the intermediate excited state in the sequential tunneling approach. We show 
that this splitting is sensitive to the number parity of electrons in the dot and the dot's ground- 
state spin. These effects are suppressed when cotunneling dominates the electrical and thermal 
conductances. We calculate the thermopower in the presence of elastic cotunneling, and show 
that some signatures of exchange correlations should still be observed with current experimental 
methods. In particular, we propose a method to determine the strength of the exchange interaction 
from measurements of the thermopower. 
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I. INTRODUCTION 

Lateral quantum dots have been studied extensively both experimentallyJ* and theoretically.— Such dots exhibit 
effects associated with both charging energy and a discrete single-particle spectrum, as captured by the constant 
interaction (CI) model. In dots with more than ~ 50 — 100 electrons, the irregular shape of the confining potential 
often leads to chaotic classical dynamics, and the fluctuations of single-particle energies and wave functions follow 
random matrix theory (RMT). The transport properties of such dots exhibit mesoscopic fluctuations as a function 
of external magnetic field and/or shape. Consequently, the electrical conductance of a weakly coupled dot exhibits 
Coulomb Blockade oscillations (a charging energy effect) and peak-height fluctuations that are well-described by 
RMT. 3-5 However, the CI plus RMT model does not explain the statistics of the observed residual fluctuations in 
peak spacings, and it was recognized that electron-electron interactions beyond charging energy play an important 
role.— A universal Hamiltonian was shown to properly describe the leading interaction terms in a chaotic quantum 
dot in the limit of large Thouless conductance^^ These interaction terms include a spin-exchange term and a Cooper 
channel term. The latter is repulsive in quantum dots and can be ignored. However, inclusion of the ferromagnetic 
exchange interaction term was shown to provide a good agreement between the calculated and observed peak spacing 
and peak height statistics in quantum dots i^r— 

The thermoelectric properties of quantum dots were studied by Beenakker and Staring 12 using the sequential 
tunneling approach in the framework of the CI model. Of particular interest is the thermopower of the dot, S = 
—AV/AT, where AV is the voltage induced by a temperature difference AT across the dot under the condition of 
vanishing electrical current. The charging energy gives rise to large periodic sawtooth oscillations of the thermopower 
(see Fig. 1) of magnitude e/2TC (where T is the temperature and C is the dot's capacitance). The thermopower 
vanishes at the charge degeneracy point where the conductance has a peak. This behavior of the thermopower 
originates in the breaking of particle-hole symmetry. For sequential transport to occur in the "valleys" between 
conductance peaks, a thermal fluctuation to an unfavorable charge state must overcome the charging gap. A valley 
of the conductance corresponds to an equilibrium state of N electrons. On one side of this valley, the N — 1 electron 
states are closer in energy (to the iV-electron states) and carry a hole current, while on the other side of the valley 
the N + 1 electron states are closer in energy and carry a particle current. In the center of the valley there is a 
(thermally rounded) discontinuity in the thermopower as it switches between hole and particle transport, leading to 
the large-scale sawtooth pattern shown in Fig. 1. The thermopower vanishes at the charge degeneracy point where 
particle and hole transport are equally likely. 

Beyond the charging energy effect on the thermopower, Ref. Il2| also predicted smaller "teeth" structures superim- 
posed on the larger scale sawtooth behavior at temperatures below the single-particle level spacing. We will refer to 
these fine structure features as quantum jumps in the thermopower. This additional structure comes from the discrete 
nature of the manifold of excited states that contribute to the thermopower. In the center of a valley, a large thermal 
fluctuation is required to transport charge through the dot. This fluctuation energy can be divided arbitrarily between 
the electron's states in the leads and states of the dot with equal probability, allowing a number of excited states in 
the dot to contribute to transport. As the Fermi energy is varied away from the valley's center, the size of the thermal 
fluctuation necessary to allow transport through the dot decreases, thus reducing the number of excited states on the 
dot contributing to charge and energy transport. The effect is that the magnitude of the thermopower increases more 
rapidly than the classical prediction between sharp negative steps that occur as each excited level falls out of the 
accessible energy range. Thus, these steps in the thermopower locate (a subset) of the excited states on the dot at fixed 
number of electrons. We note that the example discussed in Ref. [TH corresponds to an equally-spaced single-particle 
spectrum and equal level widths. Mesoscopic fluctuations have not been included in any of the theoretical studies of 
the quantum dot thermopower; these fluctuations cause the quantum jumps to deviate from uniform spacings and the 
size of these jumps to vary substantially. In Fig. lb we show a typical thermopower of the Beenakker-Staring theory 
in the presence of random matrix fluctuations in both spacings and widths, as compared with the case of an equally 
spaced single-particle spectrum and equal level widths in Fig. la. 

Experiments by Staring et alJ^- and Dzurak et alM- found evidence of the large-scale classical sawtooth structure, 
but with peak-to-peak amplitude that is significantly smaller than the theoretical prediction. In later work, Dzurak 
et alr^ found clear evidence at low temperatures of the quantum jumps predicted by the theory, but only near the 
degeneracy point and not in the center of the valleys. Subsequently, this significant discrepancy between theory and 
experiment was explained by Turek and Matveev^ who argued that cotunneling transport should dominate both 
the thermal and electrical conductances near the valley center and greatly reduce its magnitude compared to the 
predictions of the sequential tunneling theory. Furthermore, the thermopower lineshape is distorted from a pure 
sawtooth, reaching a maximum away from the valley center and then decreasing as cotunneling begins to dominate. 
Effects of sequential tunneling and cotunneling on the thermopower were also studied in single molecules that are 
coupled to metallic leads 

The Beenakker-Staring theory was formulated in the context of the CI model and did not include the effects of 
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FIG. 1. Low-temperature limit of the thermopower S (in units of S/eT) versus the effective Fermi energy If in the absence of 
exchange correlations {If is controlled by a gate voltage). The dot has 20 electrons with e 2 /2C = 55 and kT — <5/100 (S is the 
single-particle level spacing). In panel a) the level spacings and transition widths are uniform. In panel b) the level spacings 
and widths are drawn from the gaussian orthogonal ensemble of RMT. Here (and in all subsequent figures) If is taken to be 
at the degeneracy point (where the conductance has a peak). The large scale jumps at the center of the conductance valleys (at 
±5<5) correspond to a change in the excited charge state that is closest to the dot's ground state. Each sawtooth (e.g., between 
—5(5 and 5(5) in both panels contains ~ 10 fine structure (quantum) jumps corresponding to five accessible particle and five 
accessible hole single-particle levels in the dot. 

exchange correlations; hence spin entered trivially as a degeneracy of the single-particle levels. Here we are primarily 
concerned with novel effects that originate in the exchange interaction and appear in the quantum structure of the 
thermopower. The quantum structure is characteristic of sequential tunneling transport and not of cotunneling. 
The latter describes a coherent sum over many levels and smoothes out fine structure effects in the thermopower. 
The major part of this work is thus focused on the sequential tunneling thermopower, in which we find dramatic 
exchange interaction effects (in the absence of cotunneling). We then proceed to include cotunneling effects which 
impose a cut-off on the observability of these dramatic effects. This cut-off is rather stringent because the sequential 
tunneling in the conductance valleys is a thermally activated process, whereas cotunneling is not. Using current 
experimental methods and realistic device parameters, it is therefore difficult to see more than a single quantum jump 
in the vicinity of each Coulomb blockade peak. Nevertheless, we argue that the observability of any quantum jump 
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is greatly enhanced by the exchange interaction effects, suggesting that exchange might have played a crucial role in 
the observation of quantum jumps by Dzurak et al. 15 Exploiting this sensitivity to exchange correlations, we propose 
a method to determine J s from an ensemble of measured thermopower traces. 

Despite the constraints on the observability of the pure sequential tunneling thermopower, its fine structure remains 
of interest because:, unlike the conductance, thermopower directly probes the excited states of the quantum dot in 
linear response (i.e., for a small source-drain voltage). Here we generalize the Beenakker-Staring theory, formulated in 
the CI model, to a dot with electron-electron interactions beyond charging energy and in particular to the universal 
Hamiltonian framework, allowing the study of exchange interaction effects. We find that the quantum structure in 
the thermopower provides information regarding the ground-state spin of the quantum dot and the strength of the 
exchange interaction. This information can be extracted even when mesoscopic fluctuations are taken into account. 

The outline of this paper is as follows. In Section [TT] we discuss the rate equation approach for an interacting 
dot and derive an expression for the sequential thermopower in linear response theory. In particular, we obtain the 
low-temperature limit of the thermopower for the universal Hamiltonian in terms of the single-particle level widths 
and excitation energies and spins of a subset of excited states. The low-temperature results are rederived in Section 
IIIII by calculating the energy transported across the dot per charge carrier. In Section llVl we study the fine structure 
of the sequential thermopower in the presence of exchange correlations and find even-odd effects that are sensitive to 
the ground-state spin of the dot. In Section [V] we demonstrate that these effects survive the presence of mesoscopic 
fluctuations. In Section IVII we discuss elastic cotunneling in a dot described by the universal Hamiltonian and its 
effect on the thermopower. We conclude in Section IVlII with a summary and discussion of our main results. 

II. MANY-BODY RATE EQUATION CALCULATION OF THERMOPOWER 

The Beenakker-Staring model for conductance and thermopower is based on a rate equation method in which 
coherence between the dot and the leads is neglected, and the electron-electron interaction on the dot is represented 
by a constant charging energy, i.e., the CI models Hence this approach describes sequential tunneling processes. 
The linear response conductance and thermopower are given in terms of the single-particle level transition widths 
from the left and right leads and the canonical thermal occupation probabilities for the single-particle states on the 
dot. Later this rate equation approach was generalized by Alhassid et al^S, to describe transitions between arbitrary 
many-body states of the dot and applied it to the calculation the electrical conductance of a dot that is described by 
the universal Hamiltonian^ This Hamiltonian is given by 
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H = J2 £ a4 ct «a. + ^N 2 - J S S 2 , (1) 

where creates an electron in energy level e,\ with spin a, and N, S are, respectively, the total number and total 
spin operators for the electrons on the dot. The new feature in ([1]) compared to the CI Hamiltonian is the inclusion 
of the ferromagnetic exchange interaction of strength J s . Note that for an equally spaced single-particle spectrum, 
ground states of higher than minimal spin begin to occur when J s = 0.5 5, whereas full polarization of the dot occurs 
at J s — S. Mesoscopic fluctuations cause these thresholds to vary substantially from sample to sample^r— The 
effect of this exchange term has been studied in detail for the electrical conductance properties^ but not for the 
thermopower. More recently, the action for this universal Hamiltonian has been studied in a path integral approach 
and used to calculate the dot's tunneling density of states and magnetic susceptibilty . 23 i 24 

Here we calculate the thermopower in the Coulomb-blockade regime, in which the charging energy is much larger 
than the thermal energy, i.e. e 2 /C 3> kT. In this regime, we only need to consider many-body states of the dot with 
N and N + 1 electrons, where the value of N is determined by the gate voltage. We consider the case where the 
thermal energy kT is greater than the average transition width for an electron to tunnel onto the dot. This excludes 
Kondo-like resonant effects ) 25 ! 26 and allows us to use a rate-equations approach. Following Ref. [ill we consider the 
rate equations describing the time-evolution of the probability of finding the dot in each many-body states, and 
look for a steady-state solution. We expand these equations to linear order in the source-drain voltage V and the 
temperature difference AT across the dot, and arrive at a set of detailed balance equations. 

Each detailed balance equation takes the form of a sum over many-particle states; in the case of the pure electrical 
conductance, (i.e. when AT = 0), each term in the sum is separately zero^ provided that orbital occupation number 
operators commute with the Hamiltonian, as they do for the universal Hamiltonian. This term-by-term solution 
substantially simplifies the final expression for the electrical conductance. We find that for the thermopower, where 
AT =/= 0, this simplification is no longer exact, unless the transition widths are level-independent. However, when 
the temperature is much smaller than the mean level spacing (the regime of interest), we find that a term- by-term 
solution is an excellent approximation, and we use this simplification throughout this work. 
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We follow the notation of Ref. [ljj the TV electron states are indexed with i and have energies e\ N \ while the 
N + 1 electron states are indexed with j and have energies e^ N+1 \ We order them so that the ground states 
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have index 0. The equilibrium probability of finding the dot in state i is P^" ' , and similarly the probability of 



finding the dot in state j is Pj N+1 \ given by the grand-canonical statistics for the combined manifolds of TV and 
N + 1 electrons. The transition width from the TV-electron state i to the TV + 1-electron state j by an electron 

tunneling from the left (right) lead is denoted by 1^ . The effective Fermi energy in the leads, which includes 
the effect of the gate voltage, is ep, and the fractional voltage drop across the left junction is 77. We also define 
€ij = £j — £ i ~ £F an d fij = f( £ ij) where f(x) = (1 + e' 32 ') -1 is the Fermi-Dirac function with f3 — 1/kT. 
Expanding the occupancy probabilities in linear response theory as = P[ N ^ (l+^ N) /3eV + ^ N) AT/T) and 

p (N+i) = p(N+i) A + ^{N+i)p e y + $(. Ar + 1 )AT/r), we find for the thermopower: 

There is a set of detailed balance equations for the VP's and another set for the <£>'s (see Appendix A). As shown 
in Ref. 0, the equations for VP are satisfied term-by-term for a dot that is described by the universal Hamiltonian. 
This is not the case for the equations for but we find that this is an excellent approximation for T <C S. In this 
approximation, the expression for the thermopower simplifies to2£ 

s*-^^-„ — rn^. (3) 

el f 1 ij 

Z^ij - 1 i Jij r ! . .+rr. 



Using identities for the Fermi-Dirac function, and the relation Pf* = e^ eij P^ +1 , Eq. ([3]) can also be rewritten in 
the useful alternative form: 

V P N+1 ('\ - f ) r « r <j r 

S*-^^— r '-fj^- (4) 

el p N + 1 ('\ - f \ Lahi 

L-1%3 3 [ 3>r [j+ r ^ 

We next discuss the low-temperature limit. Both forms §3§ and (U]) are useful, depending on the sign of eoo, i.e., 
whether the current is particle-like or hole-like. For eoo > 0, the dot is most likely to have TV electrons, and the 
manifold of TV + 1 electron states is closer in energy than the manifold of TV — 1 electron states. Thus current will 
flow by transiently adding an electron to the dot, resulting in particle current. In this case, Eq. @ is more useful. 
At low temperatures, the largest occupation probability p^ N+1 " > of the TV + 1 electron dot corresponds to the ground 
state j = so the dominating terms in the sums of Eq. Q are the j = terms. The Fermi-Dirac function approaches 
a step function at low temperatures, and the sum over initial states i is restricted to the finite number of states with 
e i0 > 0, for which 1 — f i0 w 1 (this includes the ground state of the TV electron dot since we assume eoo > 0). Thus, 
at low temperatures and eoo > 0, the numerator of Eq. (j4]) is approximated by 

The occupation probability Pq N+1 ^ can be rather small away from the degeneracy point. However, it cancels out in 
the thermopower when expression (|S|) is combined with a similar expression for the denominator of Eq. Q to give 

1 V f , n JjoIiifl_ 
1 Z^ : e lO >o fcj 0r^+r[ n 



eT r io r ro 
z^i:e l0 >o r' +rr„ 



(6) 



For eoo < 0, the dot has TV+ 1-electrons, and the manifold of TV-electron states is closer in energy than the manifold 
of TV + 2 electron states, so the current will be hole-like. Now we use Eq. ([3]) for the thermopower to find 
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The sequential thermopower in each of expressions (J6j) and ([Jj can be relatively large, even though the electrical and 
thermal conductances are both small away from degeneracy points This thermopower depends on the gate voltage 
through the effective Fermi energy ep (which appears in ey). To see more clearly the dependence on ep, we define 

the excitation energy of the many-particle state i by e^] = — eg and rewrite eio = eoo — e ixl- ^q- © can then 
be written as 



E. (n) , e 



(N) 



Similarly, Eq. ([7]) can be written as 



eoo 1 ^<i< £ °o cx ' ir 'o+rTo , Q , 
ei ei v 1 1 io 1 m 



V (iv+i) r',.rg, 



eT eT 



■?- e ex,,- <l £ ool r D .+r 



l, V, J C J (W+l) (JV+1) , , . (iV+1) . (JV+1) W - 4.1, 4. 

where we have dehned e^ x ■ = e) e o and used eoj — eoo + e<; x j ■ Since eoo = e e — ep , we see that 

5 is a piecewise linear function of ep with a slope of 1/eT (twice the average classical slope at temperatures kT 3> S). 
However, as the effective Fermi energy varies, the number of terms contributing to the sums in Eq. ([8]) or Eq. (|9]) 
changes, and the thermopower exhibits a discontinuity or a jump. These jumps are in one-to-one correspondence with 
a subset of the many-body excited states. For gate voltages with eoo > 0, the jumps in the thermopower correspond 
to the excited states of the ./V-electron dot that have a non-zero tunneling matrix element with the N + 1-electron 
ground state. The maximal allowed excitation is e 00 (as dictated by the condition e^o > 0). For gate voltages with 
eoo < 0, the jumps in the thermopower correspond to the N + 1-electron excited states that have a non-zero hole 
tunneling matrix element with the iV-electron ground state. Such excitations are bounded from above by | eoo | (to 
satisfy the condition eoj < 0). The distance of the jump from the degeneracy point eoo — is just the excitation 
energy of the new state that appears in the sums in either (jHJ or (|9]). As |eoo| increases away from the degeneracy 
point, the number of excited states contributing to these sums increases in a stepwise manner. The maximal allowed 
value of | eoo | is ~ e 2 /2C, hence the largest excitation contributing to the sequential tunneling thermopower is half 
the charging energy. 

We focus our attention on the universal Hamiltonian, for which the orbital occupation numbers and total spin are 
good quantum numbers. In this case, the excitations that contribute to Eq. © move an electron to a level just above 
the Fermi energy and create a hole below the Fermi energy for an electron to tunnel into. Similarly, excitations that 
contribute to Eq. ([7]) move an electron from the Fermi energy to a level above the Fermi energy. The excitation 
spectrum we observe in the sequential thermopower is thus that of the single-particle energy levels (which arise from 
the non-interacting portion of the universal Hamiltonian) plus any additional spin splitting due to exchange, which 
differentiates between states with the same orbital occupation numbers. 

For a Hamiltonian that is spin-rotation invariant, such as the universal Hamiltonian, the many body states are 
characterized by good spin quantum numbers i — (aSM), where S, M are the total spin and spin projection quantum 
numbers, and a denotes all other quantum numbers. The sums over states i and j in Eqs. ([5]) and ([7]) include a 
summation over magnetic quantum numbers (including the magnetic quantum number of the ground state 0) that 
can be carried out explicitly using the Wigner-Eckart theorem to factorize out the dependence on the magnetic 
quantum numbers as a Clebsch-Gordan coefficient. For the universal Hamiltonian, orbital occupations are also good 
quantum numbers and non-zero tunneling matrix elements correspond to many-body states of N and N + 1 electrons 
that differ by the occupation of one single-particle orbital A. The tunneling width between the ./V-electron state 
i = (aSM) and the N + 1 electron state j = (a'S'M 1 ) is then given by 

T % = zshi (S M V2 ^ M ') 2 (^'||4||a5) 2 r'/ (10) 

where m = M' — M is the magnetic quantum number of the electron that tunnels into the dot, and are the widths 

of level A to decay to the left or right leads. The reduced matrix element (a'S"||ajj|aS) (which is independent of the 
magnetic quantum numbers) is given by Eq. (4) of Ref. H. Using Eq. (Ti"U|) and the unitarity of the Clebsch-Gordan 
coefficients, the sums over M, in, M' in Eqs. ([3]) and (fj| can be carried out explicitly. In particular, the particle-like 
thermopower of Eq. (J6)) is now given by 



«:ci >o ' -/ r < +rr 

Si:£ l0 >o( 2 ^ + l)r'+r{. 
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where Si is the spin of the intermediate excited state i of the ./V-electron dot characterized by an empty level A-j below 
the Fermi energy. Similarly, the hole-like thermopower of Eq. ([7]) is given by 



r' . r: 



1 Ej :£0j <o e oj( 2S 3 + 1 )rVTrf" 



S ~-^ FiTF^' (12) 

where Sj is the spin of the intermediate excited state j of the N + 1 electron dot with a single electron occupying 
level Xj above the Fermi energy. 



III. PHYSICAL ORIGIN OF THE QUANTUM STRUCTURE IN THE THERMOPOWER 

Our results for the thermopower agree with a simple physical picture we have adapted from Ref. [IH, which treats 
the thermopower in the limit of a quasi-continuous spectrum in the dot. We first use an Onsager relation connecting 
the thermopower to the Peltier coefficient II: 

„ II 1 din 

Here II is defined as the derivative of the thermal current Iq with respect to the particle current I under the condition 
of zero temperature difference AT = 0. Thus the thermopower can be determined from the heat carried by an electron 
as it is transported across the dot in a steady-state solution. 

We discuss the case eoo > 0, where the current is particle-like, eoo is the minimal energy required to add an electron 
to the dot. However, this energy can arise by any combination of thermal excitation of the electron in the leads and 
thermal excitation of the dot. Let the excitation energy of the electron in the left lead be Aj and the excitation energy 
in the dot be A 2 . These excitations will occur with probabilities proportional to e~^ Al and e~ /3A2 , respectively. Since 
the two excitations occur independently, the probability of both occurring is proportional to e ~0( A i+ A 2) . While Ai 
is, in principle, unbounded, the probability that a total thermal fluctuation is greater than eoo is negligible. Thus, 
by energy conservation, we must have Ai + A 2 = eoo and these different modes of transport are equally probable. 
However, only Ai contributes to the Peltier coefficient; this is the energy ultimately transported across the dot when 
this electron hops off to the right lead. 

Since A 2 is a dot's excitation energy, and hence quantized, the allowed values of Ai = eoo — A 2 are also quantized and 
vary with eoo- Suppose that the excitation energy A 2 corresponds to exciting the dot to a state i, i.e., A 2 = Cj ' —e$\ 
The heat carried by the electron across the dot is then Ai = e^o- To evaluate the Peltier Coefficient, we average over 
all states i that have an excitation energy less than eoo, i-e., all states i with ejo > 0. The thermopower is then given 
by 

This is exactly the expression found in the rate equation approach, under the assumption that all the transition rates 
are equal. If we now include a level- dependent weighting factor ai in the average, to allow for different tunneling rates 
into the various excited states, we find 

5= (lg) 

This gives the correct expression © for the thermopower if we identify ai = Si i0 f^° , the simplest combination of the 

partial widths rj and r^ that is both symmetric in left and right leads, and vanishes if any of the partial widths is 
zero. 

In the limit of a continuous energy spectrum and uniform transition widths, the values of Cio are uniformly dis- 
tributed in the range [0,eoo] and have an average of eoo/2. We then find S = — eoo/2eT, in agreement with the 
Beenakker-Staring result in Ref. [IH for the thermopower in the "classical" limits The case eoo < can be similarly 
treated by considering the energy required to remove an electron from the dot, i.e., the energy carried by a hole that 
tunnels onto the dot. 

At low temperatures kT « 6, the energies e^o (for a given value of ep) cannot be treated as uniformly distributed 
over the allowed range but assume values determined by the respective excitations in the dot. The condition e^o > 
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in the sum of Eq. (fT5|) is equivalent to e cx ■ < eoo- Therefore as the gate voltage (or equivalently the effective Fermi 
energy ip) is varied away from the degeneracy point, eoo increases and the number of terms in the sums of Eq. JT5j) 
increases by one each time another excited state is enclosed in the interval [0,eoo]- Each of these terms represents 
an intermediate excited state in the process of moving an electron across the dot. When a particular excited state i 
becomes energetically allowed (i.e., e^] < eoo) a $ we go further away from the degeneracy point, we observe a jump 
in the thermopower. 

IV. EXCHANGE AND NUMBER PARITY EFFECTS 

The effect of the exchange interaction is to split degenerate spin states that have the same orbital occupation 
numbers. We will show that this leads to a certain structure of the quantum jumps in the thermopower that depends 
on the number parity of electrons in the dot. This effect also depends on the ground-state spin of the dot. We 
will assume that the ground state of the odd-electron dot is S = 1/2 and discuss separately the cases where the 
ground-state spin of the even-electron dot is S = or S = 1. We note that the occurrence of an S = 3/2 ground 
state (for an odd number of electrons) is much less likely than the occurrence of an S = 1 ground state (for an even 
number of electrons) when J s < 0.5 5. 

A. Singlet ground state 

Consider an even-electron dot with a sufficiently small exchange coupling constant so that its ground state has spin 
5 = 0. When an excited state of the even-electron dot is created as an intermediate state in the tunneling process 
discussed in Sec. HH this state can be either a singlet (S = 0) or a triplet (S — 1), depending on the combined spin 
state of the two singly occupied orbitals. Since the singlet and triplet states have the same orbital occupations, they 
are split by an amount 2J S , independent of the specific single-particle spectrum. Therefore each of these states (which 
are degenerate in the absence of exchange) will appear as a quantum jump in the thermopower at values of the gate 
voltage that are separated by 2J S . In contrast, for an odd-electron dot with a ground state of spin 1/2, the allowed 
intermediate state will also be of spin 1/2 and no splitting will occur (assuming the even-electron dot has an S = 
ground state). Thus the density of the jumps on the side of even number of electrons is twice as high as the density 
on the side of odd number of electrons, with pairs of even jumps separated by 2J S . This number-parity effect is 
demonstrated in Fig. 2 for the thermopower of a dot with an exchange interaction of J s = 0.3 5 and for the case of 
equally-spaced single-particle levels and equal level widths. 

The exchange-splitting discussed above is determined by the number parity of electrons in the dot, irrespective 
of whether the process is particle-like or hole-like. The example shown in Fig. 2 describes a particle-like process on 
the even side and a hole-like process on the odd side, but similar 2 J s splitting occurs for an even-electron dot and a 
hole-like tunneling process, and there is no splitting for an odd-electron dot and a particle-like process. 

B. Triplet ground state 

The lowest S = state for an even-electron dot is described by double occupancy of its lowest N/2 levels. For a 
sufficiently small exchange coupling constant, this S — state will be the ground state of the dot. However, when 
the energy difference between the N/2 + 1 and N/2 orbitals is less than 2 J s , the dot will have an S = 1 ground state. 
In this case, the intermediate excited states in the even-electron dot can still have either S = or S = 1, leading to 
the 2 J s splitting in the quantum structure of the thermopower as discussed in Sec. IIV Al However, the intermediate 
excited states of the odd-electron dot have 3 unpaired electrons and can have either S = 1/2 or S = 3/2 spin states, 
and lead to a splitting of 3J S in the quantum structure of the thermopower. In Fig. 3 we show the appearance of 
this 3J S splitting for a dot with an equally-spaced spectrum and an exchange constant of J s = 0.6 S (for which the 
even-electron dot has an S — 1 ground state). The first jump is not split, because it corresponds to an excitation of 
an electron from the highest doubly occupied level to the singly occupied level just above, and the number of unpaired 
electrons does not change. All other jumps are split by 3J S . Thus, when the even-dot ground state is a triplet, the 
density of jumps is equal on the odd and even sides, and these jumps are paired on both sides. However, the splitting 
between paired jumps is 2 J s on the even side, and 3 J s on the odd side. 

For an equally-spaced spectrum, the even-electron S — 1 ground state occurs first (as we increase exchange constant) 
for J s = 0.5 5, and above that value there is the "Stoner staircase" of higher spin ground states, leading to even larger 
values of the exchange splitting. Typical values of the exchange constant in semi-conductor quantum dots are usually 
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FIG. 2. Thermopower for a quantum dot with equal level spacings and equal level widths in the presence of exchange interaction 
J s = 0.3 5. Results are shown for e 2 /2C = 58 at kT — 5/100. When the number of electrons in the dot is even, intermediate 
excited states can have either spin or spin 1 (see inset in upper left corner; we use the heuristic of aligned spins for spin 1 
and antialigned spins for spin 0), leading to an energy splitting of 2J a for each pair of jumps. The amplitudes of the singlet 
and triplet jumps are different because of the spin weighting factors in Eq. (JTTJ. When the number of electrons in the dot is 
odd, the intermediate excited states can only have spin 1/2, and there is no splitting. The dotted line shows the thermopower 
in the absence of exchange (i.e., J a — 0), and has been offset vertically by eT/S for clarity. For If < 0, with an even number of 
electrons the process is particle-like and the upper left inset shows excited states of the dot before the tunneling of an electron. 
For If > 0, the process is hole-like; however, for comparison with the even-electron case, the lower right inset shows an excited 
state for particle-like transport for an odd-electron dot, in which the excited states have the same spin as the ground state. 



below J s — 0.5 S. However, in the presence of mesoscopic level-spacing fluctuations, there is a finite probability to 
have higher-spin ground states at smaller values of J s , particularly S = 1 for an even-electron dot. The presence of 
3 J s splitting is a clear experimental signature of these triplet ground states. We note that the ground-state spin of a 
dot can in principle be determined by applying an in-plane magnetic field but this is a difficult experimenter— and 
level crossing at low magnetic field can lead to misidentification of the spin^ 



V. MESOSCOPIC FLUCTUATIONS 



The exchange-split quantum jumps are most easily observed at sufficiently low temperatures and when the dot's 
single-particle levels are equally spaced and have equal tunneling widths. However, in large quantum dots there are 
mesoscopic fluctuations in the level spacings and widths, and it becomes difficult to identify which jumps are paired 
together by exchange splitting, and which happen to be close to each other because of the mesoscopic fluctuations 
in the level spacings. Nevertheless, it is important to note that the splitting persist in the presence mesoscopic 
fluctuations and, given the particular spin value of the ground state, their value is independent of the particular 
sample. The jumps will be split by 2 J s on the even side of the dot, and, if the ground-state spin of the even-electron 
dot is S = 1, by 3 J s on the odd side, regardless of level fluctuations (ignoring samples for which the ground-state spin 
is larger than S = 1). Fluctuations in the transition width for each level will affect the size of the jumps, but will not 
change their position, and thus the jumps will still be separated by 2J S or 3 J,. In Fig. 4 we show the thermopower 
for a dot with the same RMT set of single-particle levels and widths as in Fig. lb, but now including an exchange 
interaction with a strength of J = 0.3 5. The greater density of jumps on the even side is apparent, but without prior 
knowledge of the energy levels, it is difficult to identify which jumps on the even side are paired. Hence a statistical 
analysis is needed. Such analysis becomes substantially more complex due to the effects of cotunneling (see Sec. IVI[) . 
and appears to be difficult to carry out with current experimental techniques. 

An interesting statistical quantity is the distribution of the spacing between neighboring jumps in the thermopower. 
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FIG. 3. Thermopower for a quantum dot with equal level spacings and equal level widths, in the presence of exchange J 3 = 0.6 S. 
All other parameters are as in Fig. 3. The 2J S splitting for an even number of electrons is still present. Now, however, the 
ground-state spin of the even-electron dot is S = 1, causing a 3J S splitting to appear for a dot with an odd number of electrons. 
The dotted line shows the thermopower for J s = 0, and has been offset vertically by eT/5 for clarity. 




FIG. 4. Thermopower for a quantum dot whose single-particle energy levels and level widths are sampled from RMT, and for 
an exchange constant of J s = 0.3 8. Results are shown for e 2 /2C = 55 at kT — 5/100. The energy levels have been labeled, 
and the inset shows the spectrum with the width of each line proportional to the conductance through that level. The dotted 
curve, offset for clarity, is the thermopower for the same RMT sample but without an exchange interaction (J 3 = 0). For this 
particular sample, the ground-state spin of the even-electron dot is a singlet, so we observe 2J S splitting on the even side and 
no splitting on the odd side. Note that the width fluctuations have removed a jump corresponding to E14. 
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FIG. 5. Histogram of jump spacings generated by sampling single-particle levels and level widths from the gaussian orthogonal 
ensemble, and evaluating the corresponding thermopower line shapes from Eq. © for an exchange constant of J s = 0.3 5. 
This distribution has two peaks at 2J S and 3J S , whose relative heights provides a measure of the probability of spin 5=1 
ground states at this value of J s ■ These peaks in the jump spacing distribution cannot be observed with current experimental 
techniques (because of cotunneling effects), but this might change if ultra- low voltage measurements become feasible in the 
future. 

These spacings will be smoothly distributed when they arise from different single-particle levels (because of level 
spacing fluctuations). However, jump spacings that arise from exchange-split jumps, are expected to lead to a large 
spike at ~ 2J S and a smaller spike at ~ 3J S . Here we assume J s < 0.5(5, so spin S = 1 ground states are not too 
frequent and higher spin ground states are rare. A histogram of this distribution is shown in Fig. 5. It is constructed 
from the thermopower line shapes of different samples drawn from the gaussian orthogonal ensemble and confirm 
our expectations. Presently, it is unrealistic to observe such a distribution experimentally, since in the presence of 
cotunneling, it is difficult to observe more than one quantum jump (see Sec. IVI|) . However, if the measurement of 
small voltages across a quantum dot can be substantially improved to allow for a more isolated dot (through which 
the conductance is smaller), then it might be possible to measure several such jumps before the contunneling cutoff 
sets in. Due to the logarithmic dependence of this cotunneling cutoff on the conductance, one would have to measure 
conductances that are several orders of magnitude smaller in order to push the cotunneling cutoff several level spacings 
away from the degeneracy point. Recent advances in using capacitively coupled quantum point contacts to measure 
extremely small currents^ suggest that perhaps such an experiment might be feasible in the future. 

VI. COTUNNELING 

The rate equation approach of Sec. |TT] assumes that the dominant transport process across the dot is sequential 
tunneling. However, at low temperatures, and away from the Coulomb-blockade conductance peaks, this is not 
the case. Cotunneling, the tunneling of electrons across the dot through virtual excitations (rather than thermal 
excitations), becomes increasingly important as the temperature is reduced. There are two cotunneling processes: 
inelastic cotunneling describes the virtual tunneling of an electron into and out of the dot that leaves the dot in a 
different state (with the same number of electrons), and elastic cotunneling that leaves the dot in the same state. 

The thermopower is given by the ratio Gt/G, where G — 9//3AV|at=o, an d Gt = dI/dAT\v=o are, respec- 
tively, the electrical conductance and thermal conductance coefficients. Making the approximation that the various 
contributions to the conductance and thermal conductance are additive, we have^i 
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In the middle of the conductance valley, £j^; c i uentlal anc [ (^sequential are p r0 p rtional to exp(— e 2 /2CkT), and are 
thus both are very small at low temperatures. If sequential tunneling is the only transport process, the exponential 
factor cancels between numerator and denominator, and the thermopower is not suppressed in the conductance valley. 
Including cotunneling increases both Gt and G. However, the relative increase in G is greater than that in Gt, and 
the thermopower is suppressed when cotunneling is significant. 

At the low temperatures kT <C 8 necessary to see the quantum structure in the thermopower, elastic cotunneling 
processes dominate inelastic cotunneling processes^ The elastic cotunneling conductance G elastlc in the presence of 
exchange correlations is calculated in Appendix B. Taking the example of an even-electron dot with an S = ground 
state, we find the average elastic cotunneling conductance to be 
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hG l G r 8 
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(17) 



where G 1 ^ = e^v^T^ is the conductance through the left (right) tunnel junction (with v& being the single-particle 



density of states in the dot per unit area), E c 



jC and we measure the effective Fermi energy relative to the 



degeneracy point. For J s — 0, Eq. (jTTJ) reduces to the known expression for the average elastic cotunneling conductance 
in the CI modeli&2£r— In Appendix B we also calculate the average elastic cotunneling thermal conductance (for a 
dot with zero ground-state spin) to be 
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In the middle of the valley ep — —E c /2, and G5i astIC = as is the case in the absence of exchange. In general, we see 
from Eqs. (JTTJ) and (TTc?|) that the effect of exchange correlations on elastic cotunneling is small. 

As we move away from the degeneracy point, cotunneling dominates the sequential tunneling, and it becomes very 
difficult to observe the quantum jumps discussed in Sec. [TXJ Thus, even though the quantum jumps are contained in 



Gt , the relevant cutoff for their observation occurs when c se q uentla i 
this condition reads (for J s = 0) 
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(19) 



For G l = G r = 10~ 3 e 2 /h, which are measurable values in current experiments, and kT — 6/15 (about the largest 
temperature at which the quantum fine structure can still be observed), we find ep ~ 0.6 8 — 0.7 8. This puts into 
question the possibility of observing pairs of jumps split by 2 J s , the simplest signature of exchange interactions. First, 
the cut-off implies that the pairs which could be observed are those nearest the degeneracy point, but, as can be seen 
from Eq. ([5]), the amplitude of the jump closest to the degeneracy point is determined by the lowest excitation energy 
in the dot and is thus quite small. Thus samples with smaller J s , which would avoid the cotunneling cutoff, will also 
show smaller jumps. Moreover such jumps will be rounded at finite T and our simulations indicate that in practice 
one cannot resolve jumps that lie within a couple of kT of the degeneracy point. While there are configurations of 
levels with which one can observe pairs of jumps, their occurrence is very rare. Thus, with current experimental 
methods, the cotunneling cutoff will make the paired jump signature of the exchange interaction difficult to measure. 

While this simplest signature is thus a challenge for future experiments, there is a less direct method for observing 
the effect of exchange correlations on the thermopower which is quite feasible with current experimental techniques. 
The net effect of the exchange interaction on the many-body spectrum is to increase the density of low-energy excited 
states. This is because higher spin-states which cost additional confinement energy are brought down near the ground- 
state by ferromagnetic exchange correlations. Alternatively the exchange interaction can make the ground state a 
higher spin state, leaving a lower spin excited state very near the ground state (as happens in the regime of the 
mesoscopic Stoner transition). The result is that the probability of observing even a single quantum jump in the 
thermopower (as has been achieved already experimentally in Ref. [iH ) is substantially enhanced by ferromagnetic 
exchange correlations. Hence it is possible to measure the probability of occurrence of observable quantum jumps in 
the thermopower and infer from that the value of J s . 

To be definite, cotunneling makes it impossible to observe jumps outside a certain cutoff C away from the degeneracy 
point (with current experimental methods C ~ 0.6(5 using kT = 5/15 and G ~ 10 _3 e 2 /H). In addition, the smallness 
and rounding of jumps near the degeneracy point put a lower cutoff c ~ 0.2 8 on how close to the degeneracy point a 
jump can be clearly seen. Thus, a jump can be observed if the actual first excited state is in the interval ~ (0.2 8, 0.6 8). 

In the absence of exchange correlations, the probability p of such an observable first excited state is p = J A pw(s)ds, 
where pw(s) = (7r /2)se^i s is Wigner's surmise for the nearest-neighbor level spacing distribution in the absence of a 
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FIG. 6. (a) The probability p of observing a quantum jump in the thermopower in the interval (c, C) = (0.2 S, 0.6 8) versus 
the exchange coupling constant J 3 /S. The curve is monotonic once the interval (shown in gray) (c/2, C/2) = (0.15,0.3 5) is 
excluded, and can thus be used to determine J s experimentally when J a is outside this interval. The exchange interaction leads 
to a strong enhancement of p as compared with its value in the absence of exchange (dotted line), (b) A histogram of quantum 
jumps that are observed within the interval (c, C) — (0.2 5, 0.6 <5) for an exchange strength of J s — 0.2 5. The large spike at 
2J a can be used to determine J s when c/2 < J 3 < C/2. (c) A thermopower trace (in the presence of elastic cotunneling) for a 
particular sample in which a triplet quantum jump (indicated by the arrow) is observed. 



b) 




magnetic field. Thus, the likelihood of observing a jump in the even-electron dot is p = e~\^. c l s )' 2 — e~T( c /' 5 ) 2 « 0.22. 
In Fig. 6(a) we show that exchange correlations enhance this probability p dramatically, by more than a factor of two 
for typical values of J s « 0.3, and as much as a factor of three for J s w 0.6, 

As already noted, the basic origin of this enhancement is the increase in the number of low-energy excited states 
due to the reduction in energy of higher spin states in the presence exchange correlations. To calculate this enhanced 
probability p for J s ^ we identify the lowest excited states of the N electron dot that have an overlap with the N + 1 
electron ground state (after an electron tunnels into the dot), and determine whether any of these excitations fall 
within the observability interval (c, C). We need to include the possibility of both singlet and triplet ground states, 
which slightly complicates the analysis. For sufficiently small J s , we can choose a random matrix spectrum and pick 
two adjacent spacings s and t (whose ensemble average is <5) to represent the first and second single-particle spacings 
above the Fermi level of the dot. If s > 2J S , then the ground state is a singlet S = 0, and the lowest excited states 
are the triplet with an excitation energy of s — 2J S and a singlet at excitation of s. In contrast, when s < 2J S , the 
ground state becomes a triplet and there are four relevant excited states with energies 2J S — s,t,2J s ,2J s + t all of 
which can cause jumps. However the 2J S excited state, which is just the singlet partner of the triplet ground state 
(i.e. same occupation numbers, but different spin), will always give a jump at a distance 2J S from the degeneracy 
point, independent of level fluctuations (as long as s < 2J S ). In fact this is a jump pair of the type we described in 
Sec. IIV1 but with the lower energy state corresponding to the ground state, making it unobservable since it carries 
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no energy through the dot. Hence for J s in the interval c/2 < J s < C/2 (i.e., 0.1 5 < J s < 0.3 S) the quantum jump 
that corresponds to the 2J B excited state is not suppressed by cotunneling effects, and a simple histogram of the 
jump location (measured from the degeneracy point) will exhibit a peak at 2 J s , allowing one to read off the value of 
the exchange constant. This is shown in the Fig. 6(b). The overall probability of observing a quantum jump in the 
interval (c, C) for this determined value of J s is given by the plot in Fig. 6(a) and can be used as a check of the result. 

In general, at larger values of J s , ground-state spins higher than S = 1 have non- negligible probability, and 
for each random matrix sample we have considered all possible states of an even-electron dot that contribute to 
the thermopower and whose excitation energy is in the interval (c, C). The probability p(J s ) (within the range 
< J, < 0.6(5) is given in Fig. 6(a). The interval 0.1 < J s < 0.3 is shown in gray; in this interval a simple histogram 
will determine J s as shown in Fig. 6(b) for J s = 0.2 5. Outside of this interval of J s , dots with ground state spin 
S = 1 will not have the 2 J s jump within the window of observability, and there will be no peak in the histogram, but 
rather a smooth distribution due to level fluctuations. However, the total probability of observing a jump can be used 
to infer the value of J s , since p is a monotonic function of J s outside of the gray-scale region, and the corresponding 
value of J s can be read off from this curve. This method for determining J s is experimentally feasible. The required 
ensemble of thermopower traces can be obtained by using finger gates and a back gate voltage. In Fig. 6(c) we show 
an example of a thermopower trace that exhibits a single quantum jump. Jumps are more likely to occur on the even 
side of the dot (this is the case in the example of Fig. 6), a feature that might be useful for determining the number 
parity of electrons on the dot with reasonable certainty by measuring the thermopower trace over several Coulomb 
blockade oscillations. 



VII. CONCLUSION 



We have shown how the presence of ferromagnetic exchange correlations modifies the fine quantum structure of the 
thermopower of a many-electron quantum dot that is described by the universal Hamiltonian. The quantum structure 
in the thermopower is sensitive to the excitation energies and transition widths of a subset of excited many-body states 
in the dot. In general, the exchange interaction splits the quantum jumps by an integer times the exchange constant, 
independent of mesoscopic level and width fluctuations. For the specific case when the ground state has spin S = for 
an even number of electrons and 5 = 1/2 for an odd number of electrons, there are twice as many jumps in the even 
valleys compared with the odd valleys, a signature of number parity. In principle, a histogram of the spacing between 
neighboring jumps can be used to measure the exchange constant in the system and also determine the probability of 
a triplet ground state. Cotunneling effects suppress these strong signatures under current experimentally realizable 
conditions, and the observation of these signatures would require the ability to measure much smaller conductances in 
almost-isolated dots. However, exchange correlations increase significantly the probability of observing any quantum 
jump in the thermopower near the degeneracy point (in the presence of cotunneling), and this is likely the reason 
that such jumps were observed. A detailed study of the distribution of the jump energy within a certain observability 
window and of the probability of occurrence of such jumps, can therefore be used to estimate the exchange constant 
J s (using current experimental methods) and determine the number parity of the electrons on the dot. Note that 
this determination does not require a magnetic field, eliminating one of the difficulties in measuring spin effects in 
quantum dots by attempting to apply a purely in-plane field. 
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Appendix A: Detailed balance equations 



The dot is described by the following rate equation; 
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In the presence of a potential difference eV and a temperature difference AT between the two leads, the Fermi-Dirac 
functions at the left and right lead are given by 
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where r\ is the fractional voltage drop across the left barrier. We are interested in a steady-state solution, i.e., 
dP.j N) /dt = dP\ N+1) /dt = for all i and j. In linear response, where V and AT are small, we expand the probabilities 



and P. 
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to first order in eV and AT around the respective (grand-canonical) equilibrium probabilities P i 
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and P, 
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p (N) = p(N)^ + ^ >f3eV + q^'AT/T) , 
p (N + i) = p(N + i )(1 + xf+Dpey + A T/T) , 

where /3 = 1/kT. Expanding the Fermi-Dirac functions to first order in eV and AT, we find 



fr. = fij -(l- ^eVf'n 



(A4a) 
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where /y — / (fcy) = [1 + e^ £iJ ] _1 , and /y denotes differentiation with respect to energy. 

Inserting Eqs. (|A4[) and (|A5I) into the right-hand sides of Eqs. (jAl|) and (|A2|) (setting the left-hand sides to zero), 
and applying further simplifications as in Ref. [l9|, we find two set of linear equations. The set for <3/'s 
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is identical to the set derived in Ref. IS The new set of equations for the <fr's is 
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The electrical current through the left lead is given by 
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and in linear response 
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Setting / = 0, we obtained our main result for the thermopower S = V/ AT 
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If the detailed balance equations (| A6|) and (| A7|) are satisfied term-by-term, the thermopower can be written as in 
Eq. ([3]). For the universal Hamiltonian, Eqs. (| A6|) are satisfied term- by-term^ but this does not generally hold for 
Eqs. dMI- 
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Another general expression for S can be obtained by calculating the current through the right lead^ For 77 = 1 
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Since the currents in the left and right leads must be equal, we can take the average of Eqs. (|A9[) (at 77 = 1) and 
(lATTj) to find 
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A simplified expression for the thermopower followed for a dot that has left- right symmetry^ i.e., T\- = r£- for all i 
and j. In this case expression (|A12|) reduces to 
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For a symmetric dot fr^r^f = 3!^ ■ , and Eq. (|A13j) is equivalent to Eq. ([3]). Thus for a dot with left- right symmetry, 
Eq. ([3]) holds generally, even when the detailed balance equations are not satisfied term-by-term. 



Appendix B: Electrical and thermal conductances for elastic cotunneling 

In this appendix we calculate the elastic cotunneling thermal conductance G§! astlc . The tunneling Hamiltonian is 

H t un = ^2 tka ^ ( r a) c M + h.c. , (Bl) 

k,a 

where V^( r a) is the dot's field operator creating an electron at the point contact r a (a — r, l) 7 t^ a is the tunneling 
amplitude for an electron in lead a with momentum k, and Ck l0 , are annihilation operators in the leads. 

Following Ref. [13, we denote the cotunneling transition rate from lead a to lead b as the dot makes a transition 
from state (N,i) to state (N,i') by W$. We have 



W$= de a r(e a ){l-f b (e b ))T$ 
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where e b — e a + e\ N ' — q. For elastic cotunneling, i = i' and e b = e a . 

We calculate the elastic transition width Tf b in second order perturbation theory 
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where E c — e 2 /C. The first contribution to the amplitude corresponds to an electron-like current; the second 
corresponds to a hole-like current. Note that — T l £. 
The cotunneling current is^ 
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where 
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Expanding the Fermi-Dirac functions f l ' r in AT and eV, we find the electrical and thermal eleatic cotunneling 
conductances to be 
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In the limit kT < 5, /'(e) ->• -5(e) and Eqs. (|B6]) reduce to 
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where the expression for G , ^J astlc is obtained after expanding T\l(e) to first order in e. The thermopower obtained from 
Eqs. (|B7D is in agreement with Mott's rule although the latter is derived for a non-interacting electron gas (see, e.g., 
in Ref. I39I). Our results hold in the presence of interactions in the dot. 

We focus on the universal Hamiltonian for which the matrix elements in Eq. (|B3p can be evaluated explicitly. As 
an example, we take a dot with even number of electrons N and a ground-state spin Si =0. The intermediate states 
j in Eq. (IB3|) for the dot with N ± 1 electrons have spin Sj = 1/2 (in the case of Si ^ 0, there are two possible values 
Sj = Si ±1/2). Since orbital occupations are good quantum numbers for the universal Hamiltonian, the sums over j 
reduce to sums over single-particle levels 
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where = 2? J2k ^( e k ~ e a)l^| 2 j 0a ( r ) is the orbital single-particle wave function A, and e\ is measured with respect 
to the Fermi energy. 

We next calculate the average of (|B8P over the mesoscopic fluctuations^. Assuming the wave function amplitudes 
at the left and right point contacts are uncorrelated and using (f>\(r a )(f>* (r a ) = S^^/A (A is the area of the dot), we 
obtain 
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where the remaining average is over the single-particle level fluctuations. Replacing the sums over single-particle levels 
by integrals we obtain 
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Using Eq. (IB10[) in expressions (IB71) . and taking « 1 away from the degeneracy point, we obtain Eqs. (|17[) 

and (|18p for the average values of the electrical and thermal conductances in elastic cotunneling. 
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